Observability of quantum phase fluctuations in cuprate superconductors 
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We study the order parameter phase fluctuation effects in cuprate superconductors near T = 0, using 
a quasi-two-dimensional d-wave BCS model. An effective phason theory is obtained which is used 
to estimate the strength of the fluctuations, the fluctuation correction to the in-plane penetration 
depth, and the pair-field susceptibility. We find that while the phase fluctuation effects are difficult 
to observe in the renormalization of the superfluid phase stiffness, they may be observed in a pair 
tunneling experiment which measures the pair-field susceptibility. 
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Underdoped cuprates exhibit large deviations from the 
predictions of BCS mean field theory, including a large 
gap A which does not scale with the transition temper- 
ature T c , and a "pseudogap" feature in the normal state 
0. These facts, together with the empirical scaling of 
the superfluid phase stiffness (SPS) and T c with doping 
p|, have motivated a picture of the cuprates in which 
the discrepancies are attributed to strong phase fluctua- 
tions H . The estimated phase fluctuation energy scale is 
smaller than the gap (pairing) energy scale in these ma- 
terials, and within this picture T c is determined by the 
small SPS rather than the gap energy. The pseudogap 
is ascribed to a precursor pairing amplitude whose phase 
coherence is destroyed above T c B. 

As many other explanations for these effects have 
been put forward, it is of great interest to devise 
tests of the phase fluctuation scenario which distin- 
guish it from others. Some evidence for thermal phase 
fluctuation effects was provided by Corson et al. ||, 
who observed unusual conductivity resonances in under- 
doped Bi 2 Sr 2 CaCu208+,5 (Bi2212) near T c and analyzed 
their data using two-dimensional Kosterlitz-Thouless- 
Berezinskii dynamics in the terahertz range, ft has been 
claimed that phase fluctuations are the dominant excita- 
tions even at low temperatures T«r c ||, determining 
the well-known linear-T dependence of the penetration 
depth 0. However, Millis et al. || have shown that 
quantum phase fluctuations cannot account for this be- 
havior. Quantum phase fluctuations do have important 
consequences for the superconductor-insulator transition 
[0], and the c-axis optical conductivity |l(]]. The above 
probes provide indirect observation of the phase fluc- 
tuation effects through electronic observables which are 
modified by the phase fluctuations. In this paper, we 
search for more direct probes of phase fluctuations within 
a model of quasi-two-dimensional cZ-wave BCS supercon- 
ductors with interlayer Josephson coupling near optimal 
doping. 

The in-plane SPS can be expressed as D a b = 
n S} abTi 2 /4md where n Sja {, is the planar superfluid elec- 



tron density, m is the effective mass of the quasiparticle, 
and d is the interplanar spacing. Because D a b is de- 
termined by the quasiparticle properties which are not 
strongly renormalized by the phase fluctuations, we find 
that although the renormalization of the Debye- Waller 
factor (e 1 *) is relatively strong, both the SPS at T = 
and its temperature corrections are weakly renormalized, 
in contrast to the case of a Josephson junction array 
(JJA) model We also consider an experiment to 

measure the excess current in a tunnel junction, which 
can be directly related to the pair-field susceptibility E^] 
X = -i0(i)([A(r,t),At(O,Q)]}. For the cuprates we find 
that this current is experimentally observable, due to the 
combination of large phase fluctuations and a low-lying 
c-axis plasmon mode. We predict a pronounced peak in 
the in the excess current at the c-axis Josephson plasma 
frequency lu c . 

We begin with a continuum BCS model with d-wave 
pairing symmetry in an isolated two-dimensional layer at 
temperature T = 1/(3 in the superconducting state: 
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with the interaction strength g > 0. Here and through- 
out the paper we set h = c = kg = 1 for con- 
venience except when numerical values are estimated. 
The gap has the property that J d 2 r e~ 4p ' r A(R, r, r) = 
A(R, r)(p 2 — pl,)/p 2 - Then we factor the pairing field as 
A(x, t) = |A(x, T)|e l ^ x ' r ); in what follows we will as- 
sume that the amplitude of the order parameter is con- 
stant, |A(x,r)| = A, and focus on the phase degree of 
freedom, 0. In order to decouple the <f> field from the 
order parameter amplitude, we perform a singular gauge 
transformation Vv( x > t) — c CT (x, T)e~ l ^' T ^ 2 , with ip a 
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the field operators for the transformed quasiparticle. The 
phase-quasiparticle coupling terms are then 



S T = d 2 x 



dr 



V0 7t i„ t , (V0) 2 ;t . . 
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where t/> is the Nambu spinor. The wavevector of the 
phase fluctuations has an upper cutoff A c of order ^q -1 ~ 
A/vf since beyond this momentum scale the mean-field 
assumption breaks down. Here we consider only the 
effect of longitudinal phase fluctuations since the pro- 
duction of vortex pairs is energetically unfavorable near 
T = and far away from the insulator transition. In or- 
der to study a realistic model, we consider such layers of 
two-dimensional superconductors with an interlayer dis- 
tance d and a weak Josephson tunneling (J) between 
adjacent layers, and a three-dimensional Coulomb inter- 
action V(q) = 4ire 2 /e b q 2 , where q = ((lu,q±), with q» 
and q± the in-plane and c-axis components of q, and 
e b the background dielectric constant. After integrating 
out the fermions we can obtain the effective phase-only 
action; the Gaussian term is 



8V(q) 



(q, w n )0(-q, -u n ) (3) 



with u> n = 2nnT the bosonic Matsubara frequencies. 
The plasma frequency to p is defined through w 2 (q) = 

{uj 2 ab q 2 +u 2 q 2 ± )/ (q 2 + q\) where u ab = ^/4irn a be 2 / e b md, 
n ab is the planar charge-carrier density at the plasma 
resonance frequency, and uj c = y/ 'Air J d 2 e 2 / e b is the c- 
axis Josephson plasma energy. Equation (^) gives the 
correct plasma spectrum for a layered superconductor, 
with u>ab the planar plasma frequency at T = which is 
0.44-1.4 eV [Q, and uj c is the c-axis plasma frequency 
which is about 0.6 meV in Bi2212 jlj] (larger in other 
materials), with an interlayer distance d « 1.5 nm and 
a dielectric constant of e b ~ 10 JlJ]. The in-plane SPS 
D a b = n s . a b/4:fnd can be read off from Eq. (|J) as the co- 
efficient of (V a &0) 2 /2 in the u) — ► and |q| — * limit. 
The quasiparticle damping term, which will broaden the 
plasma mode in Eq. (||), has been omitted, but the ef- 
fect is known to be small even in the case of the d-wave 
quasiparticle spectrum Jl4|,[]l|, and it will not affect the 
order of magnitude estimate of the correction to the SPS 
in what follows. 

We first estimate the strength of the phase fluctua- 
tion from Eq. (||). One measure of the strength is the 
Debye- Waller factor a — ', which for instance can 
be determined from the c-axis optical conductivity [fL0| . 
For our model, at T — we have 
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From this expression we see that the size of the quan- 
tum phase fluctuations is determined by the ratio of the 
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FIG. 1. The schematic diagram of the quasiparticle-phase 
field coupling and the effective action of the phase. The 
solid lines represent the quasiparticle and the dashed lines 
the phase field. 

Coulomb energy of a Cooper pair to the plasma energy; in 
the cuprates the short coherence lengths and small super- 
fluid densities conspire to enhance these fluctuations. For 
instance, assuming £o ~ 20 A and with the parameters 
given above, we estimate from Eq. (|4|) that (<p 2 ) ranges 
from 0.1 to 1, which is a sizable number compared to, 
for instance, 10 -3 in Pb. In this paper we study the ef- 
fect of these strong phase fluctuations in the BCS model 
given in Eq. ([!]), which is more appropriate near optimal 
doping, far away from the insulator transition. 

Since we are interested in the renormalization of the 
SPS, we need to go beyond the quadratic expansion in 
Eq. For instance, in the JJA model, the effec- 

tive SPS can be substantially renormalized due to the 
non-trivial potential of the form cos(0i — <f>j) fllf , In 
our model, higher-order terms can be determined by ex- 
panding Eq. (|^) and integrating out the fermions with 
a (i-wave gap. Each n-point vertex of the V^-field is a 
fermion loop (see Fig. [l]). Therefore, the renormalization 
of the SPS and its temperature dependence is determined 
by the magnitudes of the fermion loops. From the new 
effective theory of the phase fields thus obtained, 

Seff[<f>] = J2 [ dW- f d 4 X 2n T^ n \x 1 ,X2,...,X2n) (5) 

where T^ 2 ™) are the 2nth-order phason vertices, we can 
estimate the renormalization of the in- plane SPS by using 
the one loop expansion in <p as in Fig. g. At T = 0, we can 
show that the diagrams in Fig. || (a), (b) and (c) cancel 
one another in the limit of zero external momentum and 
frequency by using the identity id UJ Q(p,uj) = Q 2 (p,u>) 
and performing the integration by parts in u). Conse- 
quently, the only contribution comes from (d). The cor- 
rection to the SPS is found to be 
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Assuming that the in-plane penetration depth is 
\ ab = \J mc 2 d/4:Ttn s ^ ab e 2 ks 2000 A, we estimate that 
SD ab /D ab < 10 _1 . This should be compared to a 40% 
reduction obtained using the JJA model [ pl| . In our 
model, the phason vertices are determined by the d-wave 
quasiparticle fermion loops, which are smaller than the 
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FIG. 2. Correction to the superfluid phase stiffness in one 
loop expansion in the phase field. 

vertices of a JJA model; consequently the renormaliza- 
tion of the SPS is smaller. 

Next we study the temperature dependence of the 
SPS. The BCS theory gives a linear temperature de- 
pendence due to thermal excitations of quasiparticles, 
D ab (T) w D ab (0) - aT where a = v F ln(2)/47rv A d. Here 
t>A is the slope of the gap at the node in momentum 
space. We find that the one-loop correction to a is 



8a/ a 
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The increase in the slope a due 



to phase fluctuations is therefore hardly a measurable 
quantity, in agreement with the JJA result Ad- 
ditional temperature dependence can be obtained from 
considering classical phase fluctuations or by coupling 
the phason to a heat bath 0. However, classical phase 
fluctuation effects are not relevant to our model and the 
coupling to the heat bath leads to only a sub-leading T 2 
correction to the SPS ||. 

Unlike quasiparticle properties represented by the SPS 
as discussed above, the phase fluctuation effects on 
Cooper pair properties can be significant. Here we pro- 
pose an experiment which can measure the strength and 
the form of the phase fluctuations via the pair-field sus- 
ceptibility. We consider a c-axis tunnel junction between 
two cuprate superconductors as illustrated in Fig. || jl6) . 
The Josephson coupling between the phases (denoted 
by cj) and <fi') of the two superconductors will lead to 
the usual Josephson current oscillating at a frequency 
of 2eV/h if there is a potential difference V across the 
junction. There will also be a quasiparticle tunneling 
current. In addition, an excess current will flow due to 
the Josephson coupling of the superconducting pair-field 
of one superconducting electrode to the fluctuating pair- 
field of the other. To isolate the excess current a small 
magnetic field is applied parallel to the junction to sup- 
press the Josephson current, and the quasiparticle tun- 
neling current must be modeled and subtracted JO] . The 
excess current is interesting because it can be related to 
the pair-field susceptibility at a frequency 2eV/h Jl7| , 
and can thus provide information about the spectrum of 
phase fluctuations. 

To specialize the experiment to our case, we suppose 
that both of the electrodes are identical with a gap A 
and ignore the fluctuations in the amplitude of the order 
parameter assuming that 2eV -C A. For simplicity, we 
consider a junction in the ab plane (at z = 0) of dimen- 
sions L x x L y , with a magnetic field H y in the 6-direction, 
and we will work at zero temperature. If we assume that 



the thickness of the electrodes is larger than A c , the c-axis 
penetration depth, the Josephson coupling Hamiltonian 
for a phase difference 8<j){r, t) = 4>(r, t) — <p'(r, t) is 



d 3 r e iqxX e iS ^ r ' tS >5(z) + h.c, (7) 



L X Ly 



where Ej is the Josephson coupling energy, S 
is the junction contact area, q x « AeH y X c /hc |17| ], and 
lu = 2eV/h. The current through the junction is 



d 3 r e iqxX {e iS ^ r ' tS >)S{z). (8) 



If we calculate the current to zeroth order in Hj, and 
carry out the averages with respect to the Gaussian ac- 
tion in Eq. (|j|), we obtain the Josephson current with a 
critical current I c = (2eEj /h)a which is renormalized 
by phase fluctuations through the Debye- Waller factor a 
(there is no quasiparticle current in our model). By cal- 
culating the current to first order in Hj (linear response) , 
we obtain the excess current, 



q x ) 



eE 2 

— rjlm D R (q x ,q y = 0,w;z 
bn 



where the retarded pair field susceptibility is 
D R {r,t) = -i6(t) 



(r,t) e -i<50(O,O) 



0), (9) 



(10) 



For lu > we have Im D R (q, lu) = Im D(q, lu), where 
_D(q, lu) is the Fourier transform of the time-ordered cor- 
relation function, which for a Gaussian action is 

D(T,t) = _ ia 2 e 2(T[0(r,t)0(O,O)]> 

«- l a 2 {l + 2(T[0(r,t)0(O,O)]}}, (11) 

where the factor of two comes from the two sides of the 
junction and T is the time-ordering operator. Equation 
( |ll| ) assumes an expansion in the small parameter In a. 
If we neglect the boundary effects near the junction, we 
can obtain the propagator for the phase fields from the 
action in Eq. (|^); after an analytic continuation, we have 



(</>(q,w)0(-q, -lu)) 



-4iV(q) 



j 2 -LU 2 (q) 



i0+' 



(12) 



For an isotropic plasma frequency lu p it can be shown 
that the excess current consists of a series of (^-functions 
at integer multiples of lu p . This result is somewhat aca- 
demic since in the known isotropic superconductors the 
plasma energy is much larger than the gap energy, and 
we would expect amplitude fluctuations and quasiparti- 
cle damping to completely obliterate this effect. For an 
anisotropic plasma frequency these resonances become 
broadened — since the junction is localized at z = 0, we 
must integrate over q±, which results in a plasma fre- 
quency that ranges from lu c (when q x /q± — > 0) to Lu a b 
(when q±/q x — > 0). The result of the calculation is 
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FIG. 3. Excess pair tunneling current where uj c = 0.6 meV 
and a = 0.75. The inset figure illustrates a pair tunneling 
experiment with two cuprate superconducting electrodes SC 
and SC. We assume a junction contact area of 10 -8 m 2 , a 
normal state junction resistance Rn ~ 30 fl, and an in-plane 
magnetic field of H ~ 0.3 Gauss. [See Eq. ©.] 



4x(w,fe) 



q x e b L y L x h J ^(cj 2 _ w 2)( w 2 fc _ w 2)' 



(13) 



To quench the background Josephson current, a field may 
be chosen such that q x = 2t:/L x ; a result is shown in 
Fig. | 

The excess current exhibits a sharp onset at u> = oj c 



with a peak of the form (or 



-1/2 



This peak would 



be rounded by the small quasiparticle damping which 
we have not considered here. This experiment would 
serve as a direct observation of the phase fluctuations 
and as an alternative way to measure the c-axis Joseph- 
son plasma energy fl~i| lying below the maximum gap. 
In addition, it would provide a measure of the strength 
of the phase fluctuations. The gapless collective modes 
studied in Refs. 12 Ijj are due to the order parameter 
fluctuations which are rendered visible near the transi- 
tion temperature, which are in principle observable in any 
superconductors, whereas the quantum phason modes 
that we studied are observable as a result of the in- 
trinsically strong phase fluctuations and the quasi-two- 
dimensionality of cuprate materials. Therefore, the result 
shown in Fig. |^ is a special zero-temperature property of 
cuprate superconductors. 

We have shown that in the case of quasi-two- 
dimensional d-w&ve BCS model that we used here, 
the resulting correction to the absolute value and the 
temperature-dependence of the in-plane SPS is minute 
despite the strong phase fluctuations. Since quasiparticle 
interaction effects can obscure the fluctuation corrections 
to the SPS, it may be difficult to observe the effects of 
phase fluctuations on the penetration depth. As a more 
direct measurement, we have proposed a pair tunneling 



experiment which can probe the strength and spectrum 
of the quantum phase fluctuations. We expect that the 
pair-field susceptibility will show a pronounced peak at 
u> = lo c . It will be also interesting to explore the role of 
the order parameter fluctuations at the superconductor- 
insulator transition in the undcrdoped regime via the 
suggested experiment. In the undcrdoped regime, the 
simple BCS model fails, especially at the superconductor- 
insulator transition which is one extreme example of the 
renormalization of the SPS, and the physics of doped 
Mott insulators needs to be taken into account. 
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